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ABSTRACT 

We present the results of a iGband imaging survey of 40 arcmin 2 in fields around 
14 radio-loud AGN (6 radio galaxies and 8 quasars) with z > 1.5. The survey, 80% 
complete to K < 19.2 and complemented by i?-band imaging, aimed at investigating 
whether extremely red objects (EROs) are present in excess around high-z AGN, 
and to study the environment of z > 1.5 radio galaxies and radio-loud quasars. At 
18 < K < 19 the differential galaxy counts in our fields suggest a systematic excess 
over the general field counts. At K < 19.2, we find an excess of galaxies with R— K > 
6 compared to the general field. Consistently, we also find that the R — K colour 
distribution of all the galaxies in the AGN fields are significantly redder than the 
colour distribution of the field galaxies. On the other hand, the distribution of the 
R — K colours is undistinguishable from that of galaxies taken from literature fields 
around radio- loud quasars at 1 < z < 2. We discuss the main implications of our 
findings and we compare the possible scenarios which could explain our results. 

Key words: galaxies: active; galaxies: formation; galaxies: evolution; galaxies: pho- 
tometry; galaxies: statistics 



1 INTRODUCTION 

One of the main open questions of galaxy evolution is how 
and when the present-day massive ellipticals formed. On one 
hand, hierarchical models predict that such galaxies formed 
at moderate redshifts through merging and accretion charac- 
terized by modest star formation rates (e.g. White & Frenk 
1991). On the other hand, the so called monolithic scenario 
predicts that the present-day ellipticals formed at very high 
redshifts through a single, intense and rapid burst of star 
formation followed by a passive and pure luminosity evolu- 
tion (PLE) of the stellar population (e.g. Eggen et al. 1962; 
Larson 1975). 

From the observational point of view, the above prob- 
lems can be investigated by searching for "old", passively 
evolved spheroidal galaxies at high-z. Such galaxies would 
have the distinctive signature of very red colours due to the 
absence of star formation and to the strong K-correction 
effect. 

A population of extremely red objects (hereafter EROs) 
was indeed discovered with the combination of optical and 
near-IR imaging (e.g. Elston et al. 1988; McCarthy et al. 
1992; Hu & Ridgway 1994). EROs are found in empty sky 



* Based on observations made at the European Southern Obser- 
vatory (La Silla, Chile). 



fields (e.g. Cohen et al. 1999; Thompson et al. 1999; Yan 
et al. 2000), in the vicinity of high-z AGN (e.g. McCarthy 
et al. 1992; Hu & Ridgway 1994) and as counterparts of 
faint X-ray (e.g. Newsam et al. 1997; Lehmann et al. 1999) 
and weak radio sources (e.g. Spinrad et al. 1997). Their na- 
ture is still poorly known because their faintness hampers 
spectroscopic observations and because their colours are am- 
biguous. In fact, such colours can be due not only to an old 
stellar population, but also to strong dust reddening in a 
star- forming or active galaxy. 

A clear example of this ambiguity is provided by HR10 
(I — K ~ 6.5): its spectral energy distribution (SED) is con- 
sistent with that of an old elliptical at z m 2.4 (Hu & Ridg- 
way 1994), but optical, near-IR and submm observations 
proved that HR10 is a dusty starburst galaxy at z = 1.44 
(Graham & Dey 1996; Cimatti et al. 1998; Dey et al. 1999). 
Recent results showed that HR10 is not the only ERO de- 
tected in the submm, and that >10% of the faint submm 
population down to 850/im fluxes of a few mjy could be 
EROs (Smail et al. 1999). 

On the other hand, further observations suggested that 
also old (ages of > 2 Gyr) passively evolving galaxies pop- 
ulate the ERO class. The best examples are provided by 
LBDS 53W091 (R - K = 5.8; z = 1.55; Dunlop et al. 1996; 
Spinrad et al. 1997), ERO CL 0939+4713B (R - K = 7; 
z ~ 1.6; Soifer et al. 1999), and by most of the EROs 



© .... RAS 



2 A. Cimatti et al. 



at z ~ 1.3 around the QSO 1213-0017 (Liu et al. 2000). 
Also near-IR spectroscopy confirmed that both dusty and 
old galaxies contribute to the ERO population (Cimatti et 
al. 1999). 

Understanding the nature and deriving the abundance 
of EROs is important to shed light on the controversial issue 
of the deficit of high-z elliptical galaxies: according to some 
results, the number of galaxies with the red colours expected 
for high-z passively evolved spheroidals is lower compared to 
the predictions of passive luminosity evolution (e.g. Kauff- 
mann, Chariot & White 1996; Zepf 1997; Franceschini et al. 
1998; Barger et al. 1999). However, other works did not con- 
firm the existence of such a deficit up to z ~ 2 (e.g. Totani 
& Yoshii 1997; Benitez et al. 1999; Broadhurst & Bowens 
1999; Schade et al. 1999). 

As a first step to understand the nature of EROs, we 
started an imaging survey program aimed at selecting com- 
plete samples of such galaxies both in "empty" fields and 
around high-z radio-loud AGN. In this paper, we present 
the results of our survey around radio-loud AGN at z > 1.5. 
The main motivations of such a survey are to investigate 
whether EROs are more numerous around high-z AGN (as 
suspected in previous works; e.g. McCarthy et al. 1992; Dey 
et al. 1995), to select old passively evolving galaxies at z > 1 
by their optical/near-IR colours, to study the environment 
of high-z radio-loud AGN (e.g. Hall, Green & Cohen 1998; 
Hall & Green 1998 and references therein), and to search 
for galaxy cluster candidates at z > 1.5. We recall that the 
most distant spectroscopically confirmed cluster known to 
date has z = 1.27 (Stanford et al. 1997; Rosati et al. 1999). 
Other works selected cluster candidates around quasars at 
z > 1.3 and suggested a significant heterogeneity of the clus- 
ter galaxies, including both passively evolving old ellipticals 
as well as younger and dusty systems (e.g. Hall & Green 
1998; Liu et al. 2000). Throughout this paper we assume 
Ho — 50 kms -1 Mpc~\ fio = 1 and S1a = unless other- 
wise stated. 



2 DEFINITION OF THE SAMPLE 

The survey presented here is based on R- and T^'-band imag- 
ing. We observed totally 14 fields: 6 around radio galaxies 
(RGs) taken from the MRC sample selected at 408 MHz 
(McCarthy et al. 1997, Kapahi et al. 1998), 8 around radio- 
loud quasars (RLQs) taken from the PKS sample selected 
at 2.7 GHz (Wright & Ostrupcek 1990). We note that 
MRC1017-220 is a broad line radio galaxy (Kapahi et al. 
1998). All the targeted AGN have 1.5 < z < 2.0, with the 
exception of the quasars PKS 1351-018 (z = 3.71) and PKS 
1556-245 (z = 2.82). The radio powers at rest-frame 5 GHz 
range from 2 - 5 x 10 27 WHz" 1 for RGs to 4 - 30 x 10 27 
WHz -1 for RLQs. Table 1 lists the relevant information 
about the sample. The fields were selected according to the 
redshifts of the AGN, to their Galactic latitude (> 20°) and 
to their good observability during the telescope runs. 

We have used the predictions of the Bruzual & Chariot 
(1999) evolutionary synthesis models to define a colour selec- 
tion criterion capable to select a complete sample of EROs 
at high-z. In particular, we used instantaneous burst mod- 
els (also called simple stellar population models, SSP) with 
different redshifts of formation (zf — 2, 3, 4, 5, 6) to describe 
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Figure 1. The R—K colours predicted for different star formation 
modes, formation redshifts (zf = 2,3,4,5,6 from left to right) 
and cosmologies (assuming no dust extinction). The solid lines 
represent single instantaneous burst models (SSP). The dotted 
lines show models with star formation rate (SFR) exponentially 
decreasing with time, SFR oc exp(—t/r), with r = 0.1 Gyr and 
t = 0.3 Gyr in the top and bottom panel respectively. 

old spheroidal galaxies. In addition, we also considered mod- 
els with exponential time-scale for the star formation rate, 
SFR cx exp(-t/r), with r = 0.1,0.3 Gyr for D. = 1,0.1 re- 
spectively. Such models correspond to cases with low resid- 
ual star formation at z < 2 (i.e. < 1 M Q yr -1 for a galaxy 
with mass M ga i = 10 11 M©), and they are capable to re- 
produce the colours and the spectra of local ellipticals, as 
well as the faint galaxy optical colour and redshift distri- 
butions (Pozzetti, Bruzual & Zamorani 1996). A Salpeter 
IMF (0.1 < m < 125 Mq) and solar metallicity have been 
adopted. Models with Scalo IMF or with supersolar metallic- 
ities reach even redder colours at high-z. Optical to near-IR 
R — K colours derived from the adopted evolutionary mod- 
els are shown in Fig. 1 for two cosmologies (Ho = 50 kms -1 
Mpc -1 , D. = 0.1, 1). 

The colour selection threshold adopted in our survey is 
R — K > 6. This allows us to select old galaxies at z > 1.2 
formed at z/ > 3 in both cosmologies. It is relevant to note 
that for r >0.3 Gyr and z/ < 3, colours R — K > 6 are 
never reached. In other words, colours R — K > 6 select 
old galaxies which had a short episode of star formation in 
early cosmological epochs. For instance, assuming no dust 
extinction, a galaxy at z wl.5 with R — K > 6 should have 
z f > 3, an age >2 Gyr (Q = 1.0) or >3 Gyr (fi = 0.1) 
and SFR « 1 Moyr" 1 . 
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Table 1. The observed sample 



Object 


Type 


z 


P5Gifz 


a 


A B 


Seeing(iv) t(K) Ama(KR) Area(ftT) 


■K comp 


N{R-K> 6) 








lO^WHz" 1 




mag 


EircsGC 


SIX.* 


circinin 


fircniin 






MRC1017-220 


RG 


1.768 


0.49 


-0.95 


0.23 


1.1 


4680 


2.49 


3.08 


19.0 


3 


MRC1040-285 


KG 


1.630 


0.17 


-0.67 


0.14 


1.2 


5228 


2.88 


2.88 


19.1 


2 


MRC1048-272 


RG 


1.558 


0.44 


-0.73 


0.27 


1.3 


3780 


2.95 


3.31 


18.9 


3 


MRC1131-269 


RG 


1.711 


0.18 


-0.71 


0.24 


1.0 


4200 


2.82 


2.82 


18.9 


1 


PKS1 143-245 


RLQ 


1.950 


1.30 


-0.18 


0.22 


1.0 


3920 


2.65 


2.84 


18.9 





PKS1 148-001 


RLQ 


1.976 


2.90 


-0.48 


0.03 


1.0 


4200 


2.73 


3.00 


18.9 





PKS1148-171 


RLQ 


1.751 


0.50 


-0.30 


0.04 


1.2 


4320 


2.88 


3.49 


18.9 


0(1)* 


MRC1217-276 


RG 


1.899 


0.16 


-1.20 


0.31 


1.2 


5180 


2.75 


3.32 


18.9 





MRC1259-200 


RG 


1.580 


0.39 


-1.05 


0.24 


1.0 


4320 


3.66 


3.66 


18.9 





PKS1351-018 


RLQ 


3.707 


1.90 


+0.04 


0.11 


1.1 


4320 


2.35 


2.79 


18.9 


2 


PKS1354-152 


RLQ 


1.890 


1.80 


-0.50 


0.24 


1.0 


4200 


3.33 


3.33 


19.0 





PKS1403-085 


RLQ 


1.763 


0.63 


-0.33 


0.06 


1.0 


4200 


2.83 


2.83 


19.0 





PKS1504-164 


RLQ 


1.790 


0.40 


-0.68 


0.42 


1.0 


3744 


2.76 


2.76 


18.9 


0(2)" 


PKS1556-245 


RLQ 


2.813 


1.50 


-0.40 


0.51 


1.0 


4200 


2.90 


3.25 


19.0 






Columns: (2) RG: radio galaxy; (3) RLQ: radio-loud quasar; (4) PsGHz - radio power at rest-frame 5 GHz; (5) a: radio spectral index 
between 2.7 GHz and 5.0 GHz (5 v oc f a ) according to the NASA/IPAC Extragalactic Database (NED); (6) Ab- Galactic extinction in 
B-band; (7) t(K): total integration time in if -band; (8) Area,(KR): common area between the K- and the ij-band images; (9) Area(X): 
area covered by the ii"-band image; (10) K cornp : if-band magnitude corresponding to the 100% completeness level; (11) N(R — K > 6): 
number of objects with K < 19.2 and R - K > 6. * : the field of PKS 1148-171 contains one object with K < 19.2 and R - K > 5.75 (3cr 

lower limit). **: the field of PKS 1504-164 contains two objects with K < 19.2 and R - K > 5.35 and R- K > 5.48 (3<r lower limits). 



OBSERVATIONS, 
PHOTOMETRY 



DATA REDUCTION AND 



Near-infrared imaging was done on 1997 April 1-3 with the 
ESO/MPI 2.2m telescope with the IRAC2B camera (Moor- 
wood et al. 1992) equipped with a 256x256 HgCdTe array 
(0.506"/pixel). We used the K' filter in order to reduce the 
thermal noise. The sky conditions were photometric and the 
seeing was around 1.0". The observations were done taking a 
number of background-limited images (typically 10-14) with 
the telescope moved 10" between each image. In each tele- 
scope position, the exposures were typically 120 seconds long 
(e.g. 12 images each with an exposure time of 10 seconds). 
The data reduction was performed using the IRAF reduc- 
tion package and using the method outlined by Villani & 
di Serego Alighieri (1999). Photometric calibration was ob- 
tained observing standard stars from the Carter & Meadows 
(1995) sample. The typical night-to-night scatter in the zero- 
points was around 0.03 magnitudes. The conversion from K' 
to K magnitudes in the SAAO-Carter system was done fol- 
lowing Lidman & Storm (1995) and Lidman (1997, private 
communication) . 

Optical imaging was done in service mode from Febru- 
ary 1997 to April 1997 with the ESO NTT 3.5m tele- 
scope equipped with the SUSI camera (1024x1024 CCD, 
0.130"/pixel). The conditions were photometric and the see- 
ing ranged from 0.6" to 1.2". The observations were done 
taking 5-6 images (each with a duration of 600 seconds) with 
the telescope moved 10 arcsec between each exposure. The 
total integration time in each field was 3000 seconds, with 
the exceptions of MRC1040-285 and PKS1351-018 where it 
was 3600 seconds. The photometric calibration was obtained 
using the Landolt (1992) standard stars. The typical night- 
to-night scatter of the zero-points was around 0.05 mag. The 



correction for the atmospheric extinction was done using 
the V-band coefficients (ay) derived by the Swiss Telescope 
at La Silla during the nights of our NTT observations (see 
Burki et al. 1995), and converting them to i?-band extinc- 
tion adopting o,r = 0.55ay according to the average atmo- 
spheric extinction curve of La Silla. The data were reduced 
using IRAF. A "super" flat-field image was obtained by tak- 
ing the median of 60 images of the observed AGN fields and 
rejecting the discrepant pixels. After subtracting the bias 
and dividing for the flat-field image, the individual frames 
of each field were coadded in order to obtain the final images. 

Photometry was done using the SExtractor image anal- 
ysis package (Bertin & Arnouts 1996). The detection thresh- 
old was set to 3u of the background intensity in a con- 
tiguous area of «1 seeing FWHM. We first identified and 
extracted the objects in the A"'-band images, and then we 
cross-correlated them with the 7?-band images. The depth 
of the 7?-band images is rather homogeneous, providing a 
S/N ~ 3 for objects with R « 25 (with a 3" diameter pho- 
tometric aperture). The magnitudes were measured in 3" 
diameter apertures and converted to 6" magnitudes by ap- 
plying the aperture corrections estimated from field stars. 
Such corrections are typically in the range of 0.08-0.15 and 
0.04-0.10 magnitudes in K'- and R- bands respectively. Pho- 
tometric errors were estimated according to the poissonian 
noise from the objects and from the background. Finally, the 
magnitudes were corrected for Galactic extinction using the 
A b values derived from the extinction maps of Burnstein & 
Heiles (1982) and adopting A R = 0.56,4s and A K = 0.07^ s 
(see Table 1). 

The completeness in each K' field was estimated by 
constructing a background image for each field and adding 
random simulated objects with Moffat profiles consistent 
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with the observed seeing. The simulated objects were then 
detected with SExtractor as a function of magnitude. The 
completeness is rather homogeneous among all the fields. 
The global completeness of the whole K'- selected sample 
is 100% for K = 19.0 ± 0.1 and it decreases to about 80% 
for K = 19.2 ± 0.1. The common sky area covered by R- 
and _K"-band imaging is 40.0 arcmin 2 . At K < 19.2, only 3 
objects are undetected in _R-band (see Figure 3). 



4 if-BAND GALAXY COUNTS g 

The main uncertainty in galaxy counts comes from the star- | ? 
galaxy separation. Because of the poor pixel sampling in the z 
IRAC2B images, the morphological classifier used by SEx- 
tractor is more reliable at moderately bright magnitudes. 
Such classifier, called "stellarity index" (S), is by defini- 
tion equal to 1.0 for stars and to 0.0 for galaxies. We found 
that S > 0.9 provided a reliable classification of stellar ob- 
jects at K <17.5. We then decided to statistically subtract 
the star counts using the method outlined by Saracco et 
al. (1997): taking advantage of the good seeing and bet- 
ter sampling of the SUSI images, we extracted a subsam- 
ple of objects that were reliably classified as stellar in the 
_R-band. In order to assess the fraction of objects that are 
misclassified in the near-IR images, we computed the ratio 
R s = N s ,k /N s ,r, where N s ,k and N s ,r are the number of 
objects that SExtractor classifies as stars (i.e. with S > 0.9) 
respectively in the IRAC2B and in the SUSI images. The 
ratio R 3 was computed as a function of the if-band magni- 
tude over a range 15. 0< K <19.0. As expected, i? s ~ 1 
for 15. 0< K <17.0, indicating that SExtractor correctly 
classifies the stellar objects at moderately bright K-hanA 
magnitudes. However, the classification efficiency drops to 
R s ~ 0.5 for 17.0< K <19.0, showing that at fainter K- 
band magnitudes SExtractor underestimates the number of 
stars in IRAC2B images. 

The galaxy counts were then corrected adopting the R s 
values computed for each K-band magnitude bin. We ex- 
cluded from the counts the objects lying at < 3" from the 
edges of the images in order to avoid photometric uncer- 
tainties due to edge effects. This implied a reduction of the 
useful total area of the 14 fields from 46.4 arcmin 2 to 43.4 
arcmin 2 . 

Table 2 and Figure 2 show the galaxy counts in K- 
band as derived from the observed fields compared with the 
literature averaged counts of Hall & Green (1998) and the 
counts of the wide-field (181 arcmin 2 ) survey of Minezaki 
et al. (1998). For K < 17.5, our counts are consistent with 
the average literature counts. For 17.5 < K < 19 the counts 
seem to have a systematic excess of a factor of ~ 1.2- 1.5 
respect to average literature counts. For K > 19 the counts 
are incomplete in most of the fields (see Tab 2). In order 
to provide an estimate of the field-to-field variations of the 
counts, we do not show their poissonian uncertainties, but 
the their rms derived from the different fields. 

We also tried an alternative approach deriving the K- 
band counts by selecting objects with _R-band "stellarity in- 
dex" < 0.9 from the total area in common between the K- 
and the i?-band images (40 arcmin 2 ). Within the poissonian 
uncertainties, we obtained results consistent with the previ- 
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Figure 2. The observed differential galaxy counts (black circles). 
The error bars on our counts indicate the field-to-field rms. The 
two vertical dashed lines show respectively the 100% and 80% 
completeness limits in _R"-band. The counts at K = 19.3 are Rs70% 
incomplete. The open circles show the counts derived by selecting 
objects with il-band stellarity < 0.9 (see text). The continuum 
line represents the literature-averaged counts taken from Hall & 
Green (1998) with their rms. The triangles are the counts taken 
from the Minezaki et al. (1998) survey. 



Table 2. The _R"-band galaxy counts 



K 


Completeness 


n 


An 


16.3 


100% 


1092 


300 


16.8 


100% 


1523 


355 


17.3 


100% 


4250 


596 


17.8 


100% 


7603 


793 


18.3 


100% 


9587 


890 


18.8 


100% 


16918 


1183 


19.3 


~70% 


14907 


1112 



n: number density of galaxies deg 2 mag ; An: poissonian un- 
certainty on n. 

ous method, thus independently confirming the existence of 
an excess of counts for K > 17.5 (see Figure 2). 



5 THE MORPHOLOGIES OF EROS 

Figure 3 shows the colour magnitude diagram for all the 
objects (stars+galaxies) extracted from the 40 arcmin 2 of 
our RG and RLQ fields in common between the K- and 
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Figure 3. The colour-magnitude diagrams for all the objects 
in the radio galaxy and radio-loud quasar fields. The horizontal 
dashed lines indicates the R — K > 6 selection threshold. The 
two vertical dashed lines show respectively the 100% and 80% 
completeness limits in A"-band. The filled triangles indicate the 
objects undetected in ij-band and their 3<r lower limits in the 
R — K colours. 



the i?-band imaging. At K < 18.9 and at K < 19.2, 8 and 
11 objects with R — K > 6 are selected respectively. The 
faintest EROs have 25.5 < R < 26.0 and they are detected 
at 2a — 3a significance level in the i?-band. In addition, we 
found three objects which are clearly detected in the K- 
band at K < 19.2, but which are invisible in the i?-band. 
Their 3a lower limits on the colours are in the range of 
R - K > 5.3 - 5.8 (see Tab. 1 and Fig. 3). 

As discussed by Thompson et al. (1999) and Cuby et 
al. (1999), low mass stars and brown dwarfs can contribute 
to the population of very red objects. A visual inspection of 
the detected objects with R — K > 6 suggested a clear non- 
stellar morphology for most of them. However, in order to 
establish more quantitatively if our sample is contaminated 
by stars, we analyzed the morphology of the R — K > 6 
objects by means of their "stellarity" indices in K- and in 
7?-band, Sk and Sr. 

The stellarities in i?-band are in the range of 0.15 < 
Sr < 0.57, and clustered around Sr=0.35, indicating that 
the morphologies are rather ambiguous in the optical. This 
is not surprising because the objects with R — K > 6 are so 
faint in the optical (R ~ 25) that the morphological classi- 
fier becomes less reliable at such low flux levels. Thus, the 
morphological analysis can be better done in the A'-band. In 
the K < 18.9 sample, one object is unresolved (Sk = 0.96), 
one is ambiguous (Sk = 0.44), and six are clearly non-stellar 
(Sk < 0.15). In the K < 19.2 sample, one object is unre- 



solved, three are ambiguous (0.4 < Sk < 0.5), and seven are 
clearly non-stellar (Sk < 0.15). We noticed that the FWHM 
of the ambiguous objects is » 1.2 — 1.8 times larger than the 
seeing, suggesting that such objects may be compact galax- 
ies rather than stars. Among the three objects undetected 
in R, one is clearly non-stellar (Sk = 0.18) and the other 
two are compact (both have Sk ~ 0.85). 

In order to have more clues, we also used other two mor- 
phological parameters provided by SExtractor: the elonga- 
tion and the ellipticity, which should be respectively equal to 
1.0 and 0.0 for perfect stellar objects. In our if-band images, 
we found that, for reliable stellar objects (Sk > 0.95) their 
average values are respectively 1.10 ± 0.05 and 0.08 ± 0.04, 
whereas the objects with R — K > 6 have values in the range 
of 1.21-2.39 and 0.14-0.58 respectively. The only exception 
is the object with Sk = 0.96 which, consistently with its 
stellarity index, has an elongation of 1.15 and an ellipticity 
of 0.13. The analysis of the elongation and the ellipticity 
favours the detected objects with R — K > 6 being mostly 
non-stellar objects. 



6 THE SURFACE DENSITY OF EROS 
6.1 Our AGN fields 

The total surface densities of EROs with R—K > 6 averaged 
over the 14 observed fields derived for the K < 18.9 and the 
K < 19.2 samples are respectively 0.20±0.07 and 0.28T0.08 
objects arcmin -2 . However, we notice that there are strong 
variations of the density of objects with R — K > 6 from field 
to field, reaching values as high as ~1 arcmin -2 in some of 
the observed fields (see Tab. 1). The errors indicated in this 
section are la poissonian uncertainties. 

The average densities de- 

crease respectively to 0.18T0.07 and to 0.25T0.08 objects 
arcmin -2 if we exclude the unresolved object. Thus, the av- 
erage density of galaxies with R — K > 6 is in the range 
of 0.15-0. 18(±0. 06-0.07) objects arcmin -2 for the K < 18.9 
sample, depending on the exclusion or the inclusion of the 
object with ambiguous morphology. For the K < 19.2 sam- 
ple, the galaxy density is 0.18-0.25 (±0.07-0.08) arcmin -2 , 
but it could increase to 0.33T0.09 arcmin -2 in case the 3 ob- 
jects undetected in ii-band were galaxies with R — K > 6. 
If we exclude from our sample the densest fields of MRC 
1017-220 and MRC 1040-285, the density of objects with 
R- K > 6 and K < 19.2 decreases to 0.17 ±0.07 arcmin -2 . 



6.2 Hall & Green (1998) RLQ fields 

The average surface density of EROs derived in our fields 
has been compared to that of the Hall & Green (1998) sam- 
ple (30 fields around RLQs at 1 < z < 2 covering a total 
area of 213.7 arcmin 2 ). First of all, the r-band magnitudes 
of Hall & Green (1998) have been transformed to R using 
their relation R = r — 0.322, and their if-band magnitudes 
have been scaled to our K photometric system using the 
conversions of Carter (1993) adopting J — K — 3, consis- 
tently with the extreme J — K colors of EROs (e.g. Hu & 
Ridgway 1994). However, we note that the results do not 
change if we adopt a bluer color such as J — K = 1.5. From 
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the Hall & Green (1998) sample, we estimate that the densi- 
ties of galaxies with R— K > 6 are 0.21±0.03 arcmin -2 and 
0.29±0.04 arcmin -2 at K < 18.9 and K < 19.2 respectively. 
Such densities agree very well with those estimated in our 
AGN fields. 

6.3 Field surveys 

Our derived surface densities have been also compared to 
the results of surveys of the general field (i.e. not around 
high-z AGN) based on R- and Jf-band photometry. In this 
respect, the main reference survey is that of of Thompson et 
al. (1999) because it has the widest field coverage available 
to date (154 arcmin 2 ), and it is therefore the least affected 
by the field-to-field variations present in other surveys with 
areas about an order of magnitude smaller. Since Thompson 
et al. (1999) made use of non-standard filters (Rcadis and 
K'), we derived colour and magnitude thresholds consistent 
with our photometric system adopting the relations K — 
K' — 0.22(H — K) (Thompson, private communication; see 
also Wainscoat & Cowie 1991) and H — K = 1, and R — 
Rcadis — 0.1 (valid for B — R = 2, Huang & Thompson, 
private communication). We adopted H — K = 1 because 
the observed colors of most EROs are in the range of 0.7 < 
H — K < 1.2 (Thompson, private communication; see also 
Hu & Ridgway 1994; Spinrad et al. 1997; Stanford et al. 
1997; Soifer et al. 1999; Stiavelli et al. 1999; Yan et al. 2000). 

Excluding the 2 EROs in the Thompson et al. sam- 
ple which are spectroscopically confirmed stars, we found 
that the densities of EROs are 0.05±0.02 arcmin -2 (8 ob- 
jects) and 0.11±0.03 arcmin -2 (17 objects) at K < 18.9 and 
K < 19.2 respectively. We mention here that Daddi et al. 
(2000), who performed a survey for EROs covering a field of 
701 arcmin 2 and 447 arcmin 2 to K < 18.8 and K < 19.2 re- 
spectively, found surface densities of EROs with R — K > 6 
fully consistent with those of Thompson et al. (1999). 

We also derived the density of galaxies with R — K > 6 
from the combination of other three field surveys: McLeod et 
al. (1995) (14.7 arcmin 2 ), Cohen et al. (1999) (14.6 arcmin 2 ), 
and from the two control fields of Hall & Green (1998) (18.8 
arcmin 2 ). The if-band magnitudes of the literature samples 
have been scaled to our photometric system using the trans- 
formations of Carter (1993) and assuming an average colour 
of J — K = 3, but the results do not change if we assume 
J-K = 1.5. The r-band magnitudes of Hall & Green (1998) 
have been transformed to R as described in section 6.2. Only 
three galaxies with R — K > 6 and K < 19.2 are present 
in the joint sample of such surveys over a total area of 48.1 
arcmin 2 . The derived surface density is 0.06±0.04 arcmin -2 
which is consistent with that of Thompson et al. (1999). 

The above comparisons imply that there is a overden- 
sity of galaxies with R — K > 6 in fields around radio-loud 
AGN at z > 1.5 compared to the general field. At K < 18.9 
and K < 19.2 the average density in our AGN fields is re- 
spectively a factor of 4 and 2.5 higher than the density in the 
general field derived by Thompson et al. (1999). In particu- 
lar, if we take the field densities of Thompson et al. (1999), 
we expect 2 and 4.4 EROs with respectively K < 18.9 and 
K < 19.2 in our 40 arcmin 2 survey. The Poisson probabil- 
ity of our observing 8 and 11 EROs with K < 18.9 and 
K < 19.2 is 0.086% and 0.37%, implying an overdensity 
at 99.9% and 99.6% confidence levels respectively. However, 



it is also important to notice that strong field-to-field varia- 
tions are present in our survey, and that in some of the fields 
the density of EROs is about 10 times higher than in the 
general field (see Tab.l). Our results suggest that the strong 
field-to-field variations are due to presence of high-z clusters 
in some of the observed fields. It is relevant to mention here 
that other observations found evidence of ERO overdensities 
around high-z AGN (e.g. McCarthy et al. 1992; Dey et al. 
1995; Hall & Green 1998; Clements 2000; Liu et al. 2000). 



7 COLOR DISTRIBUTIONS 

7.1 Comparison with field galaxies 

In order to investigate the colour properties of all the galax- 
ies detected in the 14 AGN fields, we compared their R — K 
colours with those of 245 field galaxies selected from the 
samples of McLeod et al. (1995), Cohen et al. (1999), and 
Hall & Green (1998). First of all, in order to check the va- 
lidity of our field sample, we compared the average R — K 
colours of our literature sample with those of the Thomp- 
son et al. (1999) sample. For 16.8 < K < 17.8 and 17.8 < 
K < 18.8, we found average R — K colours (3.60 ± 0.11 and 
3.83 ±0.09 respectively) which are consistent with those de- 
rived by Thompson et al. (1999) in their survey (3.69 ±0.05 
and 3.83 ± 0.04 respectively). 

Figure 4 shows the colour distributions of the galaxies 
selected from the literature field sample and from our AGN 
fields. Since both distributions contain lower limits on the 
R — K colours, we used the ASURV Rev 1.2 statistical soft- 
ware package (see LaValley, Isobe & Feigelson 1992). 

As already suggested by Figure 4, we found that the 
galaxies with K < 19.2 selected in AGN fields are redder 
than those selected from the general field. In particular, all 
the statistical tests available in ASURV show consistently 
that the two distributions are different at >99.99% (> 4cr) 
significance level. Such result does not change significantly 
if we split our total sample into the RG and the RLQ sub- 
samples: the RG sample is redder than the literature field 
sample at >99.99% significance, whereas the RLQ sample is 
redder at 99.0% significance level. 

If we exclude from this analysis the four fields with the 
highest densities of EROs (i.e. MRC1017-220, MRC1040- 
285, MRC1048-272 and PKS1351-018), the result does not 
change significantly: the galaxies in the AGN fields continue 
to be redder than the field galaxies at 99.95% confidence 
level (about 3.5a). This suggests that red colour excess is 
not due to particularly overdense fields, but to the general 
population of galaxies in all the AGN fields. Such results 
strengthens the result on the overdensity of EROs with R — 
K > 6 discussed in the previous section. 

Finally, we found that the mean R — K colours esti- 
mated with the Kaplan-Meier estimator for the literature, 
(RG±RLQ), RG, and RLQ samples are respectively: 3.717, 
4.108, 4.289, and 3.921. 

7.2 Comparison with literature AGN fields 

It is also important to compare the colour distributions of 
the galaxies in our AGN field sample with other samples 
of AGN fields taken from the literature. To this purpose, 



© .... RAS, MNRAS 000, |-?? 



The density of extremely red objects around AGN 7 




R — K 
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Figure 4. The colour distributions of galaxies. Upper panel: the 
thin line represents our total sample (radio galaxies and quasars), 
whereas the thick line shows the distribution of colours of field 
galaxies. Lower panel: the thin continuum and dashed lines show 
the colour distributions of the galaxies respectively in the radio 
galaxy and in the quasar fields. Thick line: same as in the upper 
panel. 



we used the data collected by Hall & Green (1998) observ- 
ing 30 fields containing radio-loud quasars: 10 at redshifts 
1.0 < 2 < 1.4, and 20 at 1.4 < z < 2.0. The 5 GHz radio 
powers of such quasars span the range 26.0 <log P< 28.2 
WHz -1 . From the Hall & Green (1998) sample we selected 
278 and 899 galaxies with K < 19.2 respectively in the 
quasar redshift ranges 1.0 < z < 1.4 and 1.4 < z < 2.0. 
Figure 5 already shows that the colour distributions of the 
galaxies in our fields and in the Hall & Green (1998) quasar 
fields look very similar. This is confirmed by the statisti- 
cal tests which show that the colours of the galaxies in our 
(RG+RLQ), RG and RLQ fields are not significantly differ- 
ent from those of the Hall & Green (1998) fields, irrespective 
of the redshift range. The results of sections 7.1 and 7.2 do 
not significantly change if the two RLQs with z > 2 are 
excluded (PKS1351-018 and PKS1556-245). 



Figure 5. The colour distributions of galaxies. Upper panel: 
the thin line represents our total sample (radio galaxies and 
quasars); the thick continuum and dashed lines show the colour 
distributions of galaxies in fields around radio-loud quasars at 
1.0 < 2 < 1.4 and 1.4 < z < 2.0 respectively (from Hall & Green 
1998). Lower panel: the thin continuum and dashed lines show 
the colour distributions of the galaxies respectively in the radio 
galaxy and in the quasar fields. Thick lines: same as in the upper 
panel. 



distributions of the RG fields with the colour distribution 
of large Hall & Green (1998) quasar sample (30 fields at 
1 < z < 2). The statistical tests indicated that the two dis- 
tributions are not signficantly different (i.e. they are different 
at ss 1.7(7 level). The result does not significantly change if 
we add our RLQ sample to the Hall & Green (1998) sample, 
and if subsamples in different redshift bins are considered 
(e.g. 1.0 < z < 1.4 and 1.4 < z < 2.0). Again, the results 
do not change if we exclude PKS1351-018 and PKS1556- 
245 from the statistical analysis. Our conclusion is that the 
colour difference that we found between our RG and RLQ 
sub-samples is probably due to the limited statistics of our 
survey rather than to a real diversity. 



7.2.1 RGs vs. RLQs 

The results discussed in section 7.1 show that our RG 
and RLQ sub-samples have different average colours (R — 
iv=4.289 vs. fi- if =3.921 respectively). Such a difference 
is also confirmed by the comparison of their colour distri- 
butions, for which the statistical tests indicate that the two 
distributions are different at 99.9% confidence level. How- 
ever, because of the limited statistics of our survey (only 6 
RG fields vs. 8 RLQ fields), we have compared the colour 



7.3 Searching for radial trends 

Despite the rather shallow if -band depth of our sample 
and the limited sky area covered by each Jf-band image, 
we investigated the possible existence of correlations of the 
radially-averaged colour and the number density of galaxies 
with the radial projected distance from the AGN. This is 
useful to test if there is a tendency of the reddest galax- 
ies to lie closer to the AGN, and to investigate whether the 
AGN tend to be located in denser environments, as found by 
Hall & Green (1998). Figure 6 shows the radial dependences 
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Figure 6. Upper panel: the cumulative average density of galax- 
ies with K < 19.2 as a function of the radial distance from the 
AGN. The continuum line represents the cumulative number den- 
sity of field galaxies to K < 19.2 derived from the literature av- 
eraged counts of Hall & Green (1998). Lower panel: the average 
R — K colours of galaxies with K < 19.2 as a function of the 
radial distance from the AGN. The continuum and dashed lines 
represent the average R — K colour and its ±lcr dispersion as 
derived from the literature field galaxies (see section 7.1). 



of the average R — K and number density of galaxies with 
K < 19.2 in 5 arcsec wide bins. The results are strongly in- 
fluenced by the limited statistics, and no significant trends 
are visible. However, we notice that the average R — K seems 
to be higher going closer to the AGN. Should such trend be 
confirmed by deeper imaging, this would mean that the red- 
dest galaxies are preferentially located close to the AGN. 
We also note that in the bins closer to the AGN the galaxy 
number density seems to be systematically higher than the 
density of galaxies in the general field. Finally, we mention 
that our sample is not deep enough to allow us an evaluation 
of the cluster richnesses using estimators such as the Hill & 
Lilly (1991) A^o.s statistic. 



8 MAIN IMPLICATIONS 

The main finding of our survey is the discovery of an excess 
of red galaxies in fields around radio galaxies and radio- 
loud quasars at z > 1.5. Our results are in good agreement 
with the ones of the Hall & Green survey around RLQs at 
1 < z < 2. Being our survey rather shallow (K < 19) and 
based only on two-band photometry, deep multi-colour pho- 
tometry and spectroscopy are needed to unveil the redshifts 
and the nature of the excess red galaxies. With the present 



results, we can simply envisage two main possible scenarios 
to explain our finding. 

(A) The excess galaxies are at the same resdhift of 
the AGN (z — zagn) in the field, and they represent the 
"tip of the iceberg" of the clusters or the groups where the 
AGN live. This possibility is generally supported by previ- 
ous works which showed that RLQs and RGs are sometimes 
found in rich environments over a wide range of reshifts: 
0.3 < z < 1 (e.g. Yee & Green 1987; Elligson, Yee & Green 
1991; Hill & Lilly 1991; Deltorn et al. 1997), 1 < z < 3 
(Aragon-Salamanca et al. 1996; Le Fevre et al. 1996; Pas- 
carelle et al. 1996; Dickinson 1997; Yamada et al. 1998; Hall 
& Green 1998; Pentericci et al. 1998). On the other hand, it 
is relevant to recall here that radio-quiet QSOs at most red- 
shifts are found in lower density environments (e.g. Elligson 
et al. 1991; Teplitz et al. 1999 and references therein). 

A first possibility is that the observed EROs are star- 
forming galaxies strongly reddened by dust extinction (e.g. 
Cimatti et al. 1998; Dey et al. 1999; Smail et al. 1999). In 
such a case, an excess of dusty star-forming galaxies may 
suggest that high-z radio-loud AGN live in environments 
characterized by an intense star-formation activity. In this 
respect, it is relevant to mention that Ivison et al. (2000) 
recently found an excess of submm sources with possible 
ERO counterparts around the z — 3.8 radio galaxy 4C 41.17, 
suggesting an excess of dust enshrouded star formation in 
a proto-cluster environment. Such result, if confirmed and 
generalized to other cases, would be relevant to strengthen 
a scenario where high-z radio-loud AGN live in rich envi- 
ronments where strong star formation is taking place. 

Alternatively, the observed EROs could be passively 
evolving spheroidals with negligible or no ongoing star for- 
mation activity, indicating that some of the observed AGN 
may live in clusters or groups where, similarly to lower-z 
clusters, old and luminous spheroidals represent a substan- 
tial fraction of the galaxies. 

Follow-up deep imaging and spectroscopy of the most 
overdense fields (e.g. MRC 1017-220 and MRC 1040-285) 
is under way in order to verify the existence of clusters at 
z — zagn- Preliminarly, near-infrared spectroscopy of two 
of the reddest galaxies in the field of MRC 1017-220 (z = 
1.77) suggested that their redshifts are z ~ 1.5 ± 0.25 and 
that their SEDs are consistent with that of «2-3 Gyr-old 
elliptical galaxies with no dust extinction (Cimatti et al. 
1999), strengthening the hypothesis that the red galaxies 
are probably at the same redshift of the AGN. 

Assuming that the excess galaxies are ellipticals at the 
reshifts of the AGN, we estimated their rest-frame A"-band 
absolute magnitudes in case of formation redshift zf > 
zagn > 3. For the cases at 1.5 < zagn < 2.0 and adopting 
Zf — 4, we obtain that Mk are in the range of -25.5^—26.5 
and -26-^-27.5 for (SI = 1, r = 0.1 Gyr) and (SI = 0.1, 
r = 0.3 Gyr) respectively. If compared with L* K of the lo- 
cal luminosity function of elliptical galaxies (M K — —25.16, 
Marzke et al. 1998), such luminosities suggest that the 
observed EROs have rest-frame Lk ~ 1-4 — 3.4L* K and 
L K ~ 2.2 - 8. 614 for (SI = 1, r = 0.1 Gyr) and (SI = 0.1, 
r — 0.3 Gyr) respectively. 

In the case of PKS 1351-017, if the EROs in its field 
were at z — 3.7, their Mk would be in the range of -26.5-i — 
27.5 and -28^-29.5 for (SI = l,r = 0.1 Gyr) and (SI = 0.1, 
r = 0.3 Gyr) respectively, implying very high luminosities 
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(L K ~ 3.4 - 8.6L^ for f2 = 1, t = 0.1 Gyr), and probably 
too high (Lk ~ 14 — 55L* K ) in the case of (SI = 0.1, r = 0.3 
Gyr). However, we notice that if the local L* K evolves and 
increases at high-z, as suggested by some recent works (e.g. 
Glazebrook et al. 1995; Schade et al. 1999), the luminosity 
of the EROs in excess over the local L* K would be reduced. 

In this regard, it is important to emphasize that the 
presence of moderately bright galaxies with extreme red 
colours such as R — K > 6 strongly suggests the exis- 
tence a population of passively evolving "old" galaxies at 
1.5 < z < 2 whose colours can be explained only if their 
formation redshifts are zf orm > 3. Should this result be 
confirmed by the redshift measurements of the most extreme 
EROs, this would provide very stringent clues on the earliest 
epoch of massive galaxy formation (see for instance Spinrad 
et al. 1997). 

(B) The second scenario is that the excess red galaxies 
are at z ^ zagn and physically unrelated to the AGN. 

First of all, the observed EROs could be dust reddened 
starburst galaxies. For instance, HR10 provides a clear case 
of a foreground star-forming galaxy at z — 1.4 located in 
the field of a quasar at z — 3.8 (Hu & Ridgway 1994). In 
this respect, the two EROs that we detected in the field of 
PKS 1351-017 may be foreground objects unrelated to the 
quasar. 

If z > zagn > 2, extreme colours such as R — K > 6 
can be reached only if the galaxies formed at Zf > 3 for in- 
stantaneous burst models or zj > 4 for models with r = 0.1 
or 0.3 Gyr (see Fig. 1). In such a case, the detected EROs 
would represent moderately old and very luminous passively 
evolving galaxies at z > 2, and they would provide impor- 
tant clues on the existence of "aged" galaxies at high-z. Al- 
though the existence of such galaxies is in principle possible, 
the fraction of galaxies with K < 19 and z > 2 expected in 
galaxy formation models is extremely small, ranging from 
~1% in PLE scenarios (e.g. Pozzetti et al. 1996) to ~0% in 
hierarchical models (e.g. Kauffmann & Chariot 1998). This 
means that in our fields we expect at most 4 galaxies with 
K < 19 and z > 2. We notice here that at K < 19, no 
galaxies with z > 2 have been spectrospically identified yet 
in the redshift surveys of Cowie et al. (1996) and Cohen et 
al. (1999). 

The other possibility is that the red galaxies are part 
of foreground clusters or structures at z < zagn- Although 
this case is realistic, it does not easily explain why radio-loud 
AGN are preferentially seen through foreground structures 
(in order to explain the excess of EROs around radio-loud 
AGN). In this case, it is tempting to speculate whether the 
mass of such foreground structure and the geometry of the 
system can produce gravitational lensing effects on the back- 
ground AGN. A similar scenario was already suggested by 
Hammer & Le Fevre (1990). 

At a speculative level, we can envisage a scenario where 
the light coming from the AGN is gravitationally ampli- 
fied by foreground massive structures, producing a selection 
effect where the AGN that would not be selected in flux- 
limited surveys are actually selected thanks to the amplifi- 
cation of the flux due to the gravitational lensing. In this 
respect, we notice that the 408 MHz fluxes of some of the 
AGN found in the richest environments (e.g. MRC 1017- 
220, MRC 1040-285) are in the range of 1030-1090 mjy, 
i.e. only ~8-15% higher than the selection threshold of the 



MRC radio survey (0.95 Jy at 408 MHz) (see McCarthy et 
al. 1997 and references therein). This means that, assuming 
that the intrisic non-lensed flux of the sources is just be- 
low 0.95 Jy, a small gravitational amplification of the flux 
by a factor of ~10-15% would push these AGN above the 
selection threshold of the survey. The ability of clusters of 
galaxies to gravitationally lens extended sources is clearly 
demonstrated by the luminous arcs which are background 
galaxies (i.e. extended sources) distorted and amplified by 
the cluster potential (e.g. Smail et al. 1995). Although in 
our case the geometry is different, i.e. the distance of the 
lens is not much smaller than that of the source, we are only 
considering small magnification factors which could involve 
even only one radio lobe, whose size is smaller or similar to 
that of galaxies in the optical. 

To summarize, the present data do not allow us to ex- 
plore in more details the validity of the possible scenarios. 
However, according to the above discussion, we can con- 
clude that, based on their colours and Tf-band magnitudes, 
the excess red galaxies around the AGN are likely to be at 
1 < z < 2, hence favoring the scenario where radio-loud 
AGN at 1.5 < z < 2 live in rich environments, or the case 
where the red galaxies are at z < zagn- The former case is 
relevant to understand the genesis of the AGN phenomenon, 
the relationships between radio-quiet and radio-loud AGN 
and for the use of radio-loud AGN as tracers of cluster and 
massive structures of galaxies at high-z. Follow-up optical 
and near-IR spectroscopy and submillimetre photometry of 
the selected EROs is under way in order to unveil their na- 
ture and their role in the framework of galaxy formation and 
evolution. 
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